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Abstract 

The paper deals with application of GPS radio occultation (RO) measurements 
aboard CHAMP for the retrieval of tropospheric water vapor profiles. The GPS 
RO technique provides a powerful tool for atmospheric sounding which requires 
no calibration, is not affected by clouds, aerosols or precipitation, and provides an 
almost uniform global coverage. We briefly overview data processing and retrieval 
of vertical refractivity, temperature and water vapor profiles from GPS RO obser- 
vations. CHAMP RO data are available since 2001 with up to 200 high resolution 
atmospheric profiles per day. Global validation of CHAMP water vapor profiles with 
radiosonde data reveals a bias of about 0.2 g/kg and a standard deviation of less 
than 1 g/kg specific humidity in the lower troposphere. We demonstrate poten- 
tials of CHAMP RO retrievals for monitoring the mean tropospheric water vapor 
distribution on a global scale. 
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1 Introduction 



The implementation of the Global Positioning System (GP S) enabled the de- 
velopment of the rad i o occ ultation (RO) technique (e.g. iMelbourne et al 



1994; iKursinski et all Il997f ) for remote sensing of the Earth's atmosphere. 
This technique exploits atmospheric refraction of GPS signals observed aboard 
Low Earth Orbiting (LEO) satellites. Basic observable is the atmospheric ex- 
cess phase which is used for the retrieval of meteorological quantities. The po- 
tentials of GPS RO measurements for providing vertical atmospheric profiles of 
refractivity, temperature and water vapor have been demonstrated for the first 
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time by the GPS/MET experiment (e.g. Ware et al. . 19961 : Kursinski et al. 



1997). The CHAMP (Challenging Minisatel lite Payload) GPS RO experiment 



was successfully started on Feb. 11, 2001 ([Wickert et all 120011 ) and is acti 



vated continuously since mid of 2001. Considering current and planned LEO 
satellite missions (e.g. GRACE or COSMIC), GPS RO data will provide a 
valuable data base for climatological investigations and improvement of global 
weather forecasts in the future. One challenge in processing GPS radio oc- 
cultation measurements is the data analysis in the lower troposphere. The 
refractivity retrieved fr om GPS RO data shows a negative bias in relation 
to meteorological data ( Ao et al. . 2003t Beverle et al. . 20041) . which leads to 
a corresponding bias in the retrieved specific humidity. Reasons for this bias 
are GPS receiver tracking errors, uncorrected multipath in signal propagation 
and critical refraction. The application of advanced retr i eval t echniques, as 
the Full Spectrum Inversion ( FSI) method (IJensen et all 120031 ). reduces the 
refractivity bias significantly ((Wickert et al. . 120041 ) . The FSI is implemented 
to the current version (005) of the operational data analysis software at GFZ 
Potsdam, which is available since February 2004. 



2 Retrieval technique 



The retrieval of atmospheric profiles from GPS occult ation measurements has 



been described i n detail by a number of aut hors (e.g. Melbourne et all 11994 : 
Kursinski et al. . 1997 : Wickert et all I2004T ) . Briefly, the GPS measurements 
recorded by a receiver onboard a LEO (50 Hz sampling rate in case of CHAMP) 
are used together with high precision orbit information (LEO and occulting 
GPS satellite) to derive the atmospheric excess phase with millimetric accu- 
racy. These data are transformed into profiles of the ray path bending angle 
a (a), where a denotes the so-calle d impact parameter . Assuming spherical 
symmetry the Abel transform (e.g. iFieldbo et all Il971[ ) is applied to invert 
a (a) into refractivity profiles N(r), where r denotes the ray path tangential 
altitude. 



N(r) 



77.6^ + 3.73 -10 5 ^ 



The atmospheric refractivity N(r) (Eq. 1) is related to air pressure (p), air 
temperature (T) and water vapor pressure (p w ). Assuming dry air conditions, 
refractivity is direct proportional to air density and the pressure profile can 
be derived by downward integration of the refractivity profile assuming hydro- 
static equilibrium. The temperature profile is calculated consecutively by Eq. 
1. The validation of CHAMP dry temperature profiles with radiosonde data 
as well as ECMWF analyses shows a t emperature bias less than 0.5 K (RMS 



deviation 1-2 K) between 250-20 hPa (jWickert et all 12004). 
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Fig. 1. (a): Statistical comparison (zonal mean) of CHAMP refractivity (product 
version 005, 159,672 profiles, May 2001 to August 2004) with ECMWF. (b): Number 
of profiles in relation to profile availability at 10 km altitude corresponding to left 
panel. 



However, the dry air assumption is not valid over wide areas of the mid and 
lower troposphere, leading to an ambiguity of the dry and wet refractivity 
term (Eq. 1). To deal with this problem, additional meteorological informa- 
tion is necessary. There are diffe r ent re trieval techniques using such informa- 
tion. [Gorbunov^n^^ofelovski^ (1993) describe an iterative algorithm which 
uses external temperature information (e.g. from ECMWF analyses) to sepa- 
rate dry and wet part of the measured refractivity. In result, pressure and 
humidity profile are derived from the refractivity data. Another approach 
to retrieve both humi dity and temperature profile is the lDvar technique 
(He alv and This optimal estimation method requires background 

information (temperature, humidity and pressure) as well as error character- 
istics of the measurement (refractivity) and the background (e.g. ECMWF). 
A further approach, in the following referred to as direct water vapor pressure 
(DWVP) retrieval, has been developed at GFZ Potsdam. Here background 
(ECMWF) temperature and pressure information are used to calculate water 
vapor pressure (p w ) directly form refractivity data applying Eq. 1. The dif- 
ference between the so derived p w and background humidity information is 
used to adapt the background pressure for a recalculation of p w . The pressure 
values converge very quickly and the procedure is stopped after the second 
iteration step. The known negative refractivity bias in the lower troposphere 
(Fig. 1(a)) states a general problem for the humidity derivation. To avoid 
DWVP retrieval outliers especially in that region, deviations between back- 
ground and retrieval humidity are restricted to the double ECMWF humidity 
error of the current profile during the iteration. 



At GFZ Potsdam both, lDvar and DWVP algorithms are implemented for 
tropospheric water vapor retrieval. In the following we present results from 
both techniques and compare these with radiosonde measurements. 
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Fig. 2. Comparison of vertical specific humidity (a), temperature (b) and refractiv- 
ity(c) profiles derived from CHAMP DWVP and lDvar retrieval with radiosonde 
Lindenberg and ECMWF data. Example for occultation 226, October 25, 2002, 
22:27 UTC, 51.44°N, 14.66°E. 



3 Results and validation 



To give an impression on the retrieval results, Fig. 2(a) shows an example of 
CHAMP lDvar and DWVP specific humidity in comparison to radiosonde and 
ECMWF data. Both CHAMP retrievals come to quite similar results reveal- 
ing significant improvement of the background (ECMWF) specific humidity in 
comparison to radiosonde data. Nevertheless, lDvar shows a slightly smoother 
result than DWVP. This has been observed in several cases and can be con- 
sidered as a general difference between both retrievals. Obviously, DWVP is 
more sensitive to vertical structures in the input refractivity than lDvar. Fig. 
2(b) reveals good agreement between radiosonde and ECMWF temperature 
profiles. The lDvar temperature shows only small deviations from the back- 
ground. Finally, radiosonde, ECMWF and CHAMP refractivity profiles are 
given in Fig. 2(c). Especially between 600 and 800 hPa the CHAMP refrac- 
tivity shows better agreement to the radiosonde observation than ECMWF. 
This obviously corresponds to improvements of the CHAMP humidity retrieval 
above the background humidity in this altitude range. 

The lDvar and DWVP retrieval results have been validated with radiosonde 
data. Fig. 3 (a)-(c) show the statistical comparison (bias and standard devia- 
tion) of vertical specific humidity profiles from ECMWF, lDvar and DWVP 
with coinciding radiosonde profiles on a global scale. For the years 2002 and 
2003 about 13,400 coincidences have been found (see Fig. 3 (e), coincidence 
radius: 300 km spatial and 3 hours temporal). Radiosonde data were quality 
checked by comparison with ECMWF and have been ignored in case of more 
than 10% refractivity deviation. As can be seen from Fig. 3(b) and (c), the 
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Fig. 3. Statistical comparison (years 2002-2003) of vertical specific humidity profiles 
from global radiosonde stations with: (a) ECMWF, (b) lDvar, (c) DWVP. Corre- 
sponding comparison of CHAMP and radiosonde refractivity is shown in (d). Blue 
lines represent bias, red lines standard deviation. Number of compared data points 
(coincidence radius: 300 km, 3 hours) is shown in (e). 

statistical comparison of lDvar and DWVP with radiosonde data comes to 
quite similar results. The negative CHAMP refractivity bias (Fig. 3(d)) leads 
to a specific humidity bias of of about -0.2 g/kg in the lower troposphere. 
Nevertheless, the standard deviation is similar to the result from radiosonde 
comparison with ECMWF (Fig. 3(a)). By comparing Fig. 3(a)-(d) it has to 
be noticed that between 600 and 500 hPa DWVP shows a slightly stronger 
relation to the refractivity data than lDvar which seems to be more influenced 
by the background (ECMWF) data in this region. 



4 Global application 



The CHAMP humidity profiles may be used for investigation of mean seasonal 
or medium term water vapor distribution on a global scale. Due to rather low 
data exploitation in the lower troposphere (see Fig. 1(b) and Fig. 3(e)) and 
accuracy restrictions at low humidity levels in the upper troposphere, the mid 
troposphere region is most appropriate for such investigations. Fig. 4(a) shows 
the mean global water vapor distribution at 500 hPa pressure level derived 
from DWVP results for the northern summer season of 2002 according to 
a grid of 2.5° resolution in latitude and 5.0° in longitude respectively. The 
corresponding data coverage (profiles per pixel) is shown in 4(b). It has to 
be mentioned that current and planned LEO RO missions like GRACE and 
COSMIC will provide a significantly extended data base which will allow for 
global coverage within much shorter time scales than the CHAMP mission 
alone. Furthermore, RO data will be of growing interest for climatological 
investigations in the future if first medium and long term RO data sets become 
available. 
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Fig. 4. (a): CHAMP DWVP mean global water vapor distribution at 500 hPa for 
northern summer conditions (June-August 2002). (b): Data coverage according to 
(a) (12575 data points). 

5 Conclusions 

The lDvar and DWVP techniques state valuable tools for the humidity re- 
trieval from occultation refractivity measurements in the mid and lower tro- 
posphere. Statistically, both methods come to comparable results. Even if 
a unique separation of dry and wet refractivity components is not possible, 
the derived humidity information may be improved above the background 
(e.g. ECMWF) provided that temperature background and refractivity mea- 
surement are of sufficient quality. Advanced retrieval techniques like the FSI 
method significantly reduced the negative refractivity bias in the lower tro- 
posphere. Nevertheless, the refractivity retrieval needs further improvement. 
Potentials of CHAMP RO data for global water vapor monitoring have been 
demonstrated. The CHAMP RO experiment generates the first long-term RO 
data set. Following satellite missions like GRACE and COSMIC will signif- 
icantly extend the RO database improving the capabilities for global water 
vapor monitoring and climatological investigations. 
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